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1 INTRODUCTION 



Almost every physical theory, electromagnetism, gravity and etc., is treated as the the- 
ory of a constrained system. The theoretical framework of the constrained system was 
invented in 1950's by Dirac with the Hamiltonian formalism [00. The essence of the 
Dirac method is as follows. Due to the Hamiltonian formalism, we normally use the 
canonical Poisson brackets in the unconstrained phase space. In the constrained phase 
space, however, a new type of brackets - Dirac brackets - must be introduced a la Dirac. 
In passing from the classical to the quantum theory, we replace these Dirac brackets by 
the commutators (xl/ih). To obtain the Schrodinger equation, we have to find the rep- 
resentation of the commutation relations, i.e., we must obtain the representation of the 
operator variables. 

The classical theories of constrained systems are completely described in Hamiltonian 
formalism with the Dirac method^]. However, when we pass to the quantum theory, we 
meet some problems 00]. One of such problems is that Falck and Hirshfeld presented 
a simple example of the Dirac method applied to the dynamics of a point particle They 
quantized a system with a particle on a sphere in the 3-dimensional Euclidean space. In 
this case the Dirac brackets of the variables, expressed in the Cartesian coordinates, are 
not canonical. Here "• • • canonical" means that the Dirac bracket between a coordinate 
variable and its conjugate momentum variable is 1 and that other Dirac brackets with 
respect to the variables vanish. As is clear from this example, the Dirac brackets among 
the conjugate variables are generally not canonical, and the difficulties arise to find out 
the representation of the operator variables when we pass to the quantum theory Q. 

Homma et al. suggested a new method to overcome these difficulties for a constrained 
particle system on a sphere and on a general hypersurface 00. A constraint f(q % ) = 
with q\ the coordinates of a particle, is evidently equivalent to df/dt = with a vanishing 
constraint at some instant. They defined the system under the constraint f(q l ) = 
as system(I) and that under the constraint df/dt = as system(II). The two systems 
should be equivalent both in the classical mechanics and in the quantum mechanics. 
The Dirac brackets in system(I) are not canonical, but those in system(II) are canonical, 
whichever is the reason why they introduced system(II). They, therefore, easily found 
the representation of the operators in the latter system. Now let us call these constraint 
conditions introduced initially to make a particle constrained on a specific manifold, by 
the name ^initial constraint conditions" . 

Homma et al. dealt with the system constrained only with a single initial constraint 
condition, i.e., the system constrained on a hypersurface. In other words, they formulated 
the particle dynamics in the constrained coordinate space whose dimension is only reduced 
by one from that of the unconstrained coordinate space: in the constrained phase space 
the dimension is reduced by two from that of the unconstrained phase space. However, 
the system with several initial constraint conditions has been discarded up to now. 

Following their method, we, in this paper, treat several initial constraint conditions to 
complete the dynamics of a particle on various manifolds. Moreover these conditions are 
made to depend not only on the coordinates but also on the velocities. Thus if we succeed 

1 Its details are explored in Appendix 
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to construct the theory considering the several initial constraint conditions depending on 
coordinate and velocity variables, we are to have the general theory of constrained systems 
with arbitrary initial constraint conditions. 

This paper is composed as follows. In Sec. |2|, we shall review the method of Homma et 
ai, restricted to an example of a particle constrained on a sphere of arbitrary dimensions. 
In Sec. 0, we shall deal with several initial constraint conditions depending only on the 
coordinates, and in Sec. (|, with these depending both on the coordinates and the velocities 
of a particle. The last section will be devoted to the summary and discussions. 

In present paper, the summation convention, according to which summation over the 
repeated indices is meant, is adopted. 



2 QUANTIZATION ON A SPHERE 

Now, we consider, as system(I), a nonrelativistic particle (mass : m) constrained on an 
(N — l)-dimensional sphere (radius : \fA and iV > 2) in the iV-dimensional Euclidean 
space. The initial constraint condition in this system is f(q) = q l qi — A = 0, and the 
Lagrangian is 

U = lmq% + X(q l q t -A) (z = 1, 2, • • • , N) , (1) 

where 

q ~ dt ' 

the q % denote the coordinates of the particle and A is a Lagrange multiplier. The conjugate 
momenta are defined as 

Pi = dqT = mq > ^ 
dLi , . 

OA 

Equation (|3]) is a primary constraint condition 

0i = Px ~ , (4) 

where " ~ " represents Dirac's weak equality. Following the Dirac method we obtain the 
secondary constraint conditions 

2 = q'qi-AnO, (5) 

03 = q l Pi^0, (6) 

04 ee p'pi + 2m\q t q i w . (7) 



By Eq.(0), the variable A is written in terms of q % and pf. 

P l Pi 
2mq t q i 



p l Pi 
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We, further, eliminate p\ and A from this system with Eqs.(f|) and (H), respectively. There- 
fore the variables of this system are q l and Pi, and the remaining constraint conditions 
are (j>2 ~ and 03 ~ 0. The Dirac brackets of the variables are calculated as follows: 



{q\q J }n 

{Pi,Pj}o 



0, 

A 



(9) 
(10) 

fill 



We find that the Dirac brackets (P)— (pT|) are not canonical. After we have introduced the 
Dirac brackets, the extended Hamiltonian is 



Hie = —p l Pi 
2m 



(12) 



In the canonical quantization method, we replace the Dirac brackets (|9j)- (|llf) by the 
commutators (xl/ih) and regard the variables as operators: 



q\q> 



1 ,Pj 



Pu Pi 



a Is) - m 



ih 



PiQj ~ PjOi 
A 



(13) 
(14) 

(15) 



where the symbol with caret " A " designates an operator. 

In order to obtain the Schrodinger equation, we must find the coordinate representation 
of pj. In our simple case, it is possible to obtain its representation by the trial-and-error 
method. However, since the commutation relations of q l and pi are complicated in general 
cases, it is difficult to do that. Therefore we discuss the quantization method a la Homma 
et al. 

Firstly, we consider the constrained system with the initial constraint condition q l qi = 
and refer it to system(II). This condition is the time derivative of the initial constraint 
condition in system(I): f(q) = 2q l <ji = 0. The Lagrangian of system(II) is given by 



The conjugate momenta are 



IT E 

n A E 

The primary constraint condition is 



dq 1 
dL u 



mqi + Xqi , 
. 



(16) 

(17) 
(18) 



Xi 



n, 



(19) 
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and the secondary constraint condition is 



X 2 = q% - Xq r q t « . (20) 
In system(II), the Dirac brackets of q l and IT are 

(21) 
(22) 
(23) 

Equations (PT|)-(|2"3|) show that the Dirac brackets are canonical in system(II). The ex- 
tended Hamiltonian is 



{Mn 


= o ; 


{q\^h 


= 


{n f , n 3 -} D 


= . 



1 

2m 1 



(F^Y , (24) 



where 



F\ = (25) 
q 2 

q 2 = q l qi ■ 

From Eq. (p5|) , we find that F\ is the projection matrix which projects an arbitrary vector 
onto the sphere. 

In the above discussion, we treat (q\ IT, A, IIa; i = 1,2, N) as the variables of sys- 
tem(II). However we can reduce them to (q l ,Hi;i = 1,2, ...,N) by solving the equations 
Xi ~ and \2 ~ to represent IIa and A in the other variables. In this case, since all 
constraint conditions are disappeared, we can deal with the system through the Poisson 
brackets constructed with {q l , IT; % = 1,2,..., N). The Poisson brackets of q l and IT take 
the same form as in Eqs.fl2"T|)-(|2"3"|). 

Secondly, we introduce new momentum variables P{. 

Pi = F{Uj. (26) 

which are interpreted as projected momentum variables. The Dirac brackets including Pi 
are calculated as follows: 

{g\P,} D = 5)-^f = F), (27) 
{Pi,Pjh = Piqj ~ PjQi . (28) 



Investigate the time development of q , 



da 2 

-^ = {q 2 ,H UE } D = 0, (29) 
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and we find that q 2 is a constant of motion. Putting q 2 = A at t = as the initial 
condition of the differential equation (|29|), the Dirac brackets fl27|) and ( f28|) take the same 
form as in Eqs.(^) and (|TT1). We further obtain, from Eq. (|26|) , the constraint condition 
whose form is the same as that of Eq. (13) : 



q'Pi = q l F\^ 3 = . 



(30) 



Finally, we regard Pi as the momenta Pi of system(I). In this way the constraint struc- 
ture of system(II) is equivalent to that of system(I). Therefore system(I) and system(II) 
have the same content in the classical theory. 

In order to quantize system(II), we replace the Dirac brackets (|2l"l)-(^3l) by the com- 
mutators (xl/ih) of the operators: 

(31) 
(32) 
(33) 

From these commutation relations, we immediately have the coordinate representation of 
EL: 



q\ q> _ 


= o, 


g\ft/ 


= ihS) 


EL, ft/ 


= . 



EL 



-%h 



d_ 

dq { 



(34) 



Therefore it is easy to obtain the Schrodinger equation of system(II) in the coordinate 
representation. 

Requiring the equivalence between system(I) and system(II) in the quantum theory 
and making use of Eq . (|34]) , we obtain the Schrodinger equation for system(I). 

To sum up, in the quantum case, we adopt q 2 = A as the initial condition, and then 
regard Pi as being equal to pi after symmetrizing the order of operators in Eq . (p6|) : 



Pi 



I, + II,/-', 



(35) 



We, of course, have all constraint conditions of system(I) based on those of system(II) in 
quantum level. Substitution of Eq.([3^) into the quantum Hamiltonian of system(I) 



1 

2m 



PPi , 



(36) 



gives us the Hamiltonian of the system(II). The resultant quantum Hamiltonian of sys- 
tem(II) is reduced to, with an appropriate ordering, 



8m v 



(37) 



System(II) thus becomes equivalent to system(I) also in quantum level. In order to obtain 
the coordinate representation of the Hamilton operator of system(I), we substitute Eq. 
into Eq.(]37|). 
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3 SEVERAL CONSTRAINT CONDITIONS DEPEND- 
ING ONLY ON COORDINATE VARIABLES 



In this section, we treat the dynamics of a nonrelativistic particle of mass m in an N- 
dimensional Riemannian space with a Riemannian metric Qij{q). Let a particle be con- 
strained with M initial constraint conditions 

/,(</) = (s = l,2,...,M;M<JV), (38) 

where s is the index to distinguish various initial constraint conditions. Namely, the 
particle is constrained in an (N — M)-dimensional subspace. Here we assume that the 
constraint conditions depend only on the coordinate variables q % and that these conditions 
are differentiable with respect to q l . The Lagrangian L\ in this system is expressed as 
follows: 

Li = ^m/,//V + AM)q l ~ V(q) + X s f s (q) = 1, 2, • ■ ■ , N) , (39) 

where 

9 ~ dt ' 

and V(q), Ai(q) and X s are the scalar potential, the vector potentials and the Lagrange 
multipliers, respectively. The momenta conjugate to q l and A s are 

Pi = -q^ = mgijtf + Ai(q) , (40) 
Ps = ^ = 0. (41) 
The primary constraint conditions are 



<Pla 

We thus have the total Hamiltonian 



= P' s « . (42) 



2m 

where 



Hit = ^-9 tJ (p, ~ Ai){ V] - A 3 ) + V(q) - X s f s (q) + v% , (43) 



v s = X s 



and V s are the multipliers for Eq. (f42|) . For the consistency that a constraint condition 
should hold also after the time development, we succeedingly have secondary constraint 
conditions 



i 2 , = /,(?) «0, (44) 
m oq l 



>:;, = -^^(Pj-A^^O, (45) 
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Ms 



d_ 

dq k 



d 



dq l ) 3 dq k 



(Pi - At) 



2d? 9 dqi PkPl + d/ 9 Al)9 dq* Pk 2 9 dq* dq^ 9 MM 



1 nij 9fsdV 1 
9 a i a ■> H 







(46) 



where 



D 



st = 9 



>j 



df s df t 



dq i dqi 

Since the initial constraint conditions f s = are independent, i.e., the hypersurfaces 
f s = are independent each other, the normals to the surfaces df s (q)/dq l are linearly- 
independent. This means that the determinant of {D st ) does not vanish, and X s are 
rewritten as the functions of q l and Pi with Eg . (f46|) . Moreover v s are determined as the 
functions of q l and pi for the consistency of 4s . Hence all the physical quantities are 
expressed with q l and p^. We calculate the Poisson brackets among <p as (a = 1, 2, 3, 4): 



{(p2s, 4>3t} 



— D 

m 



St ) 



-D 



m 



st J 



and all other Poisson brackets vanish. Hence all <; 
conditions. The Dirac brackets among q % and pi are 



(47) 
(48) 

are second class constraint 



0. 

df 



i\ st dfs_ dft 
dq k dqi 



-l\st 



_d_ 

dq j 



ki 



dq k 



d_ 

dqi 



v dq k . 



(49) 
(50) 

(<<-j')»(51) 



{q\q J }D 
{q\Pj}d 

where | /=o denotes that we impose all initial constraint conditions f s = after having 
calculated the Poisson brackets. We exchange the indices % and j in the first term and 
the term thus made is designated by (i <-> j). With the Dirac brackets, the extended 
Hamiltonian is 

H m = ^~g ij (pi - Ai)(jPj - Af) + V(q) . (52) 
2m 



f=o 



We now consider system(II) following the procedure in Sec. 
Lagrangian Lu is 



In system(II), the 



Lu = -mg l3 q^ + A t (q)q l - V(q) + X s f s (q) (z, j = 1, 2, • • • , N) . (53) 



The momenta conjugate to q % and A s are 

dL u 



n, 
n' 



dq* 
dLu 
dX s 



mg^q 3 +A^ + X 



dq { 



. 



(54) 
(55) 
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Therefore the primary constraint conditions are 

xis = K « o . 

Thus the total Hamiltonian is 



Hi 



IIT 



2m 



g" [ik-Ai-X 



s df± 
dq i 



"> 1 



+ V(q)+u s W s 



(56) 



(57) 



where 



u s = X s 



and u s are the multipliers for Eg . ([56"D . The consistency requires the following secondary 
constraint conditions: 



- 1 H d f* ( tt a \t 8 ^ 







(58) 



Eliminating X s and w s in the same way as in system(I), all the physical quantities are 
expressed with q % and Ilj. We calculate the Poisson brackets among Xas (a = 1,2): 



{Xu,Xn} = —D st , 
m 



(59) 



and all other Poisson brackets vanish, indicating that all Xas = are second class con- 
straint conditions. The Dirac brackets among q % and Ilj are 



{g\H,} D = 5), 

{n^n^o = o. 



The extended Hamiltonian is 



Hue = —g i] F k z (U k - A k )F l j(ILi - A,) + V(q) 
2m 



where 



F k i = 6 h i -g u [D- L r^^ 



,-i\ st dfs_ dft 
dq l dq 1 



Here we define new momentum variables Pf 

P - A = F k t (U k - A k ) . 
The Dirac brackets including Pj are calculated as follows: 

{q\P J } D = F), 



{Pi,Pj} D 



dfs 
dq 1 



(D 



-l\st 



d 



g 



M 9 ft 



dqi V dq 



Pi 



d 



M 9 ft 



dqi V dq 



(60) 
(61) 
(62) 



(63) 
(64) 

(65) 

(66) 
(67) 



9 



Moreover we have the time development of f s 



= {fs, Hue}d = 



(68) 



which shows that f s are the constants of motion. When we choose f s = as the initial 
conditions of the differential equations (|68|) (i.e., at t = 0), the Dirac brackets fl66|) 
and fl67|) take the same forms as in Eqs. (|50|) and (|5T|) . Note that since (Pi — Ai) are 
the tangential components of (U k — A k ) to the surfaces and df s /dq l are the normal 
components, (Pi — Ai) just satisfy the constraint conditions (^) with pi, replaced by 
Pi. The argument mentioned above means that the constraint structure of system(I) is 
equivalent to that of system(II). We take f s = as the initial conditions and identify Pi 
with the momenta pi in system(I). 

We pass to the quantum theory under the requirement of the above equivalence of both 
systems(I) and (II). For the hermiticity, we need to symmetrize the order of the products 
of operators expressing observables. Firstly, we symmetrize the momentum operators: 



Pi-A i = l[F k i (f = 0),IL k -A k ] 



(69) 



where [ , ] + denotes the anticommutator. The commutation relations including IT are 

~q\ flj] = ihd) , (70) 
"IT, ft;] = 0- ( 71 ) 

These are canonical. The coordinate representation of ftj is 



n,- = — ih 



Here we use the Christoffel symbol 



d_ 

J dq i 



-ih 



k 

ki 



(?)■ 



(72) 



J 

ki 



(?) = -j9 dn (dk9in + dig kn - d n g u 



(73) 



Secondly, we obtain the quantum Hamiltonian by symmetrization. We write down the 
Hamiltonian, based on Eqs.( 



and (|69|) , as follows 



Hm = ^~[F k i(f = 0),II fc - A k } + ^[F)(f = 0),IL, - A t } + + V(q) + ti 2 Q 



where 



2m 



Q = — 

2m 



F\F l j [ g 



i.i 



a 

a k 



-,3b 



I 

kb 



+\F k A { y 



ij 



a 
a I 



b 

bk 



- 0.F) 
d k F\\ g v{F\ 



d l + F\g ij (d k F k 3 ) di 



a 
a I 



+ diF l j 



(74) 



(75) 
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The quantum mechanical potential Q, typical of quantum mechanics, is a term which 
makes the Schrodinger equation be covariant under general coordinate transformations. 
The hermiticity of the Hamiltonian ([74]) leads to the vanishing of the antisymmetric parts 
for those multiplied by the metric g % K Hence we have the hermiticity conditions 

-F h i d k F l j + (d i F h j )F l h + F k j d k F\-{d j F h i )F l h = 0, (76) 

diF h jd k F k h - d j F h i d h F k h = 0, (77) 

(5 k i - F k i}d k {(S l j — F l j)Ai(q)} - (z <-> j) = 0, (78) 

F^idkF)) = 0, (79) 

Note that Eq.(|79|) comes from the hermiticity of the quantum mechanical potential. With 
these conditions, we finally obtain the Schrodinger equation 

ih— = H^l 
at 

= -^F\{F l ^j). k + l ^F\A k g^F l ^ x + ^-{F k i A k g i ^ F^) 

+ ^-{F\A k gVF l 3 A l )m + V(q)V , (80) 

where ^ is a wave function, and .k denotes the covariant derivative with respect to q k , 
and its concrete expression for T\ is 

Ti k = V k Th = d k T\ + j ^ 1 T\ -\ l ik \ Th ■ (81) 



4 SEVERAL CONSTRAINT CONDITIONS DEPEND- 
ING ON COORDINATES AND VELOCITIES 

In this section we consider the system of a particle constrained with several initial con- 
straint conditions, depending not only on the coordinates q l but also on the velocities 
g*(z = l,2,---,iV). 

We proceed by taking an example, namely by putting concrete forms for the initial 
constraints f s (q,q)- Other forms for f s (q,q) are, of course, allowed. For instance, f s (q,q) 
may include the terms whose powers of velocities are greater than one: <y' si j k ...(q)q l q-'q k • • •■ 
But this case leads to a complicated nonlinear equation of motion. We thus start with a 
simple (and general) example: 

fs(q, q) = a sl (q)q l + (3 s (q) = (s = 1, 2, • • • , M; M < N) , (82) 

where both a si (q) and fl s (q) are the functions of the coordinates q\ The conditions ( |32"D are 
assumed to be independent with respect to the velocity variables. We thus immediately 
have det(E st ) ^ with E st = g^a si a t j. (See the discussion leading to det(-D^) ^ in the 
previous section.) 
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The Lagrangian of this system is given by 

L = ['nii./l'Y + MqW ~ V (<l) + A'/.Gz, q) (s = 1, 2, M) , (83) 

with X s , the Lagrange multipliers. The momenta Pi conjugate to q % and Il s conjugate to 
X s are 

Pi = mg^q 3 + A^q) + X s a si , (84) 
IT = 0. (85) 

From Eq . (|85|) , the primary constraints are deduced: 

Xls = IT^0 (s = l,2,...,M). (86) 

Therefore the total Hamiltonian of the system is given by 

H T = ^9 ij (Pi 'A- X s a sl ) (pj - A, - X l a t3 ) - X s (3 S + V{q) + v s U s , (87) 

where v s (= X s ) are the Lagrange multipliers for the primary constraints (|5£|). 

The consistency for the primary constraints leads us to the secondary constraints 

X 2s = --E st X l + -g ij (pi - A { )a sj + (3 S « . (88) 
m m 

We can write A s in terms of q l and pi as follows: 

A s « (E-YWKPi - A)a tj + m(3 t ] . (89) 

The multipliers v s are determined as the functions of q l and Pi from the consistency for 
\2s- Hereafter we eliminate Ii s and A s with Eqs.(86) and (|89|) , and regard q l and pi as 
the independent variables of the system. The Poisson brackets of xis and X2s are 

{Xis,Xn} = -E s t , (90) 
m 

which shows that all the constraints are second class. 
Now the Dirac brackets for the variables become 

(91) 
(92) 
(93) 

The total Hamiltonian ( |STD is reduced to the extended Hamiltonian 

H E = ^9 ij \Pi -A- (E-Ya*{™Pt + 9 M (Pk - A k )a tl }} 

x [ Pj - Aj - (E- l ) st a sj {ml3t + g k \p k - A k )a a }] 

- (3 s (E~ 1 ) st [m(3 t + g kl (p k - A k )a tl ] + V(q) . (94) 



{q\q J h 


= 0. 


{q\Pj}v 


= 5 { - 

3 


{Pi,Pj}v 


= . 
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The classical dynamics of the system is described with this Hamiltonian and the Dirac 
brackets 

In order to quantize the system, we must replace the Dirac brackets (|9T1) -(|93"1) by the 
commutators (xl/ih) and regard variables q l and Pi as operator variables q l and pi. The 
results are 

(95) 

(96) 
(97) 

The commutation relations (p5|)-(p7D are canonical, so that we easily obtain the coordi- 
nate representation of the momentum operators. Symmetrizing the order of operators in 
Eq.(P^D to guarantee the hermiticity of the quantized Hamiltonian, we have 



q\ q> 


= o, 


q\Pj_ 


= ihS) 


\pi,Pj 


= 0. 



where F\ is defined by 



\F\ , p k - A k \^\F\, p,- M\ + 

8m 

- ^F\{E- V ra u fi v , p k - A k } + 



— 1\UV 



2 

- mf] u (3 v (E 
+ h 2 Q + V(q) 



F\=5\-g^ k {E 



-l\st 



f=0 



(98) 
(99) 



and Q is the quantum mechanical potential mentioned in Sec. [| Henceforth we express 
all the operators in the coordinate representation and omit the caret . 

The requirement of hermiticity further leads us to the following two conditions: 



F\F l j <f 



a 

ak 



kb 



{E- 1 ) uv a U] (3 v g i \5 k l + F k i ) 



I 

Ik 



. 



(100) 
(101) 



The Schrodinger equation of this system is written in the coordinate representation of 
the operators, with the conditions ( |100| ) and ( |101| ), as 



1 



2m 



g lj F k t (-ihV k - A k )F l {~ihV l - A t )V 



gV{E- l )™a uj (3 v F%-ihV k - A k )m 
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+ Uh[V k (E- l ra uj (3 v g^{5\ + F k 8 p 

+ V(q)V , (102) 

^ : a wave function 
where Vfc is the covariant derivative with respect to q k , defined in Eg. ([81]) . 

5 SUMMARY AND DISCUSSIONS 

We have, in this paper, presented a general theory of a particle constrained with several 
initial constraint conditions, thus aiming at completion of the dynamics of a particle on 
a general manifold. 

First we put the initial constraint conditions as functions of q l . The Dirac brackets 
of q l and pi are complicated in structure, which leads to some difficulties associated with 
quantization. Therefore we apply the method of Homma et al. to the system and obtain 
the Schrodinger equation. In this process, requiring the Hamiltonian to be hermitian, we 
have the conditions (ff6|)-([79|) . This is the result that was not obtained by Homma et al. 
and is due to the existence of several initial constraints in the system. 

Next we treated the system of a particle constrained with the initial constraint con- 
ditions (|S2T), depending not only on the coordinates but also on the velocities. The Dirac 
brackets of variables simply become canonical, so that we have the coordinate represen- 
tation of the momentum operators straightforwardly. We have the conditions (|100|) and 
( |101| ) by requiring that the Hamiltonian should be hermitian. 

Since F\ depends on f s , we regard the hermiticity conditions as the conditions for 
the initial constraints f s . Hence, in the quantum theory , the structure of the manifold on 
which the particle exists is to be restricted. 

Note that there is a case that the initial constraint conditions depend explicitly on 
time. In this case we cannot express the time evolution of physical quantities only with 
the Dirac bracket. We must make, at the same time, an explicit use of the Poisson bracket, 
therefore we cannot apply the method of the present paper to this system. The formalism 
of the quantization of this system is a subject for future study. 
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